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Abstract 

We calculate the screening charge density distribution due to a point charge, such as that of a 
positive muon placed between the planes of a highly anisotropic layered metal. In underdoped 
hole cuprates the screening charge converts the charge density in the metallic-plane unit cells in 
the vicinity of the (i + to nearly its value in the insulating state. The current-loop ordered state 
observed by polarized neutron diffraction then vanishes in such cells, and also in nearby cells 
over a distance of order the intrinsic correlation length of the loop-ordered state. This in turn 
strongly suppresses the loop-current field at the fi + site. We estimate this suppressed field in 
underdoped YE^CusOe+z and La2_ x Sr z Cu04, and find consistency with the observed 0.2-0.3 G 
field in the former case and the observed upper bound of ~0.2 G in the latter case. This resolves 
the controversy between the neutron diffraction and juSR experiments. The screening calculation 
also has relevance for the effect of other charge impurities in the cuprates, such as the dopants 
themselves. 
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In metals composed of metallic planes with negligible hopping from one plane to the 
next, the screening length for a point charge placed between the planes depends on the 
distance c between the planes, the distance d of the charge to the planes, and the electronic 
compressibility of the planes. Through an elementary calculation based on the Thomas- 
Fermi method, we obtain the screening charge distribution, and discuss how it is modified 
in strongly correlated metals. 

The immediate occasion for doing such a calculation is that sensitive muon spin relaxation 
(/iSR) experiments in underdoped YBa 2 Cu30 6+:c and La2_ x Sr x Cu04 [4j detected at 

most only a small fraction of the magnetic field expected at the muon site given the magnetic 
order observed bypolarized neutron scattering j^, [f], Q], and in the case of La 2 _ a; Sr. r Cu04 
gave a null result Q. Such a field is expected from proposed loop currents (si, [l| as illustrated 
in Fig. [H 
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FIG. 1: The orbital current pattern in the Cu02 plane. 



The discrepancy is huge: the expected field at the muon site is of order 100 G, while the 
/iSR experiments yielded 0.2-0.3 G in YBa 2 Cu306+a; [3J] and no observable field > 0.2 G in 
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arized neutron diffraction exper- 
5] have recently been verified in 



La 2 _ a; Sr :! .Cu04 

Loop-current magnetic order has been observed by po 
iments; the initial results on underdoped YBa 2 Cu 3 6+:E 

an independent experiment [a], and magnetic order of the same symmetry is also observed 
by neutron diffraction in underdoped HgBa2Cu04 [7[. In both of these compounds the ob- 
served order is long-range, limited only by the resolution of the experiment. In underdoped 
La 2 _ x Sr x Cu04 only short-range order with a correlation length of about 10 A is observed 
(around the same Bragg spots and with similar integrated intensity as the other compounds) 



10] , due presumably to the larger charge disorder in this family of cuprates and the resultant 



effects discussed in this paper 



11] . The magnitude of the ordered moment is large, estimated 



to be 0.05-0.1 hb per O-Cu-0 plaquette at about 10% doping, decreasing to zero at the 
quantum critical point; the sub lattice magnetization per unit cell is twice this value. Earlier 
dichroic ARPES results 1^] on underdoped Bi s Sr 2 CaCu 2 8+< 5 samples revealed exactly the 
same symmetry of magnetic order. These experimental results were considered controversial 
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15], but the neutron diffraction experiments lend them further credence. In under- 



doped YBa 2 Cu30 6+;c , a small dichroic signal characteristic of a time-reversal breaking order 
parameter has also been observed [16] and may be associated with a small ferromagnetic 
order accompanying the main loop-current order l?j]. 

On the basis of all these experiments, there are grounds now to believe that the loop- 
current order is universal in underdoped cuprates. It is therefore important to understand 
why /xSR experiments see little or no sign of the accompanying magnetic order. 

The loop-current order can occur only in a limited range of electron density, from about 
5% to about 15% deviation of charge density from the half-filled band; below ~5% doping 
the materials are antiferromagnetic insulators, while the quantum critical point at which the 
loop current order ends is at ~15% doping. The experiments are done with positively charged 
muons which unlike polarized neutrons are not innocent observers of a magnetic field. 
Their charge must be screened in the metal by repelling holes from their vicinity with an 
integrated density equal to that of the fi + charge. We show that this makes cells in the 
metallic planes in the vicinity of the muon acquire a local charge density characteristic of 
the insulator. This of course removes the loop-current order locally. In addition to this 
effect, the magnitude of the loop-current order parameter is reduced around such cells over 
a distance of the order of the intrinsic correlation length of the ordered state (several lattice 
constants). We calculate the field expected at the fi + site when these effects are included, and 
find that it can be of the order of the observed values in YB^CusOe+x [l, 2] and below the 
observed upper limit in La2- x Sr x Cu04 4]. Our elementary calculation and conclusion also 
has applications to other properties, such as local one-particle spectra measured by STM and 
the effect of charged defects on transport properties in layered metals and superconductors. 
We note that such screening does not occur in an insulator, where there is a gap to low 
energy charge excitations. Thus antiferromagnetism in the insulating state of the cuprates 
is perfectly well observed in /zSR experiments. 
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We use the model of metallic planes specified by the charge density of a typical under- 
doped cuprate, with each plane having a compressibility dn/dfi and a separation c between 
planes as shown in Fig. [2J 




FIG. 2: Schematic of Cu02 planes and the fi + site. 



The positive charge is placed at a distance d < c above one of the planes. We use the 
Thomas-Fermi method to calculate the potentials and the charge density changes. The 
potential <p(r, z) approximately follows the Poisson equation 

- V 2 0(r,z) = 4vr-(p ext (r,2;) + p ind (r, z)), (1) 

where e has been introduced to account for static polarizability in between the planes. The 
long wavelength limit of e is known to be about 4. The external charge is 

p cxt (r,z) = 5(z-d)5(r), (2) 

and the induced charge is given in the Thomas-Fermi approximation by 

e 2 

p ind (r,z) = — j(dn/dfi)(/)(r,z). (3) 
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It is straightforward to show that the solution of these equations is such that the charge 
density p n (r) on the nth plane can be obtained by numerically evaluating 

" w/c dk „ 



Pn{r) 



J ^TT J-ir/c ^ 



t/c 

where the radial distance r is measured from directly below the muon, 

/, e x 2K 2D cr(k,£) 



1 k 2D sinh£c 



£ cosh£c — cos kc 
1 e~ lkc smh£d + sinh£(c — d) 



2£ cosh £c — cos kc 
and the response k 2 d is defined by 

/t 2D = 2n(e 2 /e)(dn/dp). (4) 

In the general case, the screening charge distribution p n {r) can only be obtained by 
numerical evaluation. Some features of interest may however be noted immediately For a 
charged impurity near an isolated two dimensional conducting plane, the screening charge 
density decays as 1/r 3 [18]. For large enough k 2 d (of the order or larger than the effective 
mass enhancement divided by e of about 1), the screening charge density is independent of 
k 2 d and is determined solely by c and d. For the usual situation where d is close to c/2, 
the screening is determined entirely by c. For k 2 d 0, the incompressible limit as in an 
insulator, there is no screening. 

The calculated p n {f) i n the planes in the immediate vicinity of the muon is shown in 

Fig. m 

More distant planes have negligible screening charge. A check on these calculations is 
that the integrated screening charge equals in magnitude the p + charge. We have taken 
/t2D ~ 10/c, estimated from an effective mass ~ e ~ 4. The results are insensitive to the 
value of k 2 d beyond about 10/c. For much smaller values of k 2 d than this, the screening 
distance goes up with a corresponding decrease of charge density immediately above and 
below the muon. This is as it should be, since for k 2D — > we should approach the insulating 
case. 

There is an interesting general point about screening in strongly correlated material 
that is not directly included in the Fermi-Thomas approximation. In a parameter range 
as discussed below, the Thomas-Fermi screening charge in the immediate vicinity of the 
charged defect is such that the charge density crosses over past the charge density in the 




FIG. 3: The screening charge density in the planes in the Thomas-Fermi approximation for a muon 
centered between planes (d = c/2). (a): Nearest-neighbor plane, (b): Next-nearest-neighbor plane. 

insulator to effective local electron doping. However, the compressibility of the insulator is 
zero, so that the screening charge in a cell cannot go beyond the insulating regime into the 
electron doping regime. Instead, the charge density will locally stay close to zero doping, 
and the screening charge will be distributed more evenly outside the cell in the immediate 
vicinity as well as increased in more distant planes. Within the above scheme this may be 
simulated by a local response K2T>[p(r)], which goes to zero when p(r) is in the insulating 
regime. This has the effect of flattening the charge distributions compared to those shown 
in Fig. [3] and therefore increasing the size of the screening charge. 

We now apply these calculations to the parameters of underdoped La 2 _ a; Sr :r Cu04. In 
this compound the planes are separated by c ~ 6.7 A and the in-plane lattice constant is 
a »s 3.8 A. From Fig. [3j for the case d = 0.5c in the middle of the two planes), the 
integrated screening charge density in the cell immediately above and below the fi + changes 
by about 10%, and in each of the neighboring 8 cells it changes by about 5%. There is 
negligible screening charge farther away. At a typical underdoping concentration of, say, 
10%, this places the cells directly underneath definitely in the insulating regime and their 
neighbors close enough that the proposed orbital order is strongly affected in them from 
this effect alone. As noted above, incompressibility of the insulating regime increases the 
size of the screening region beyond this estimate. In addition, the suppressed loop-current 
order requires a significant additional distance to recover as shown below. Thus even the 
conversion of one cell to the insulating composition is enough to make a large /z + field 
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unlikely. 

For d = 0.3c, which is a large asymmetry in the placement of the muon, the screening 
charge density in the nearest and the next nearest neighbor plane are shown in Fig. HI 
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FIG. 4: The screening charge density in nearby planes in the Thomas-Fermi approximation, with 
a large asymmetry in the placement of the muon {d = 0.3c). (a): Nearest neighbor plane, (b): 
Next nearest neighbor plane. 

In the Thomas-Fermi approximation the charge in the nearest cell in the closest plane 
changes by about 30%, and the charge in in the 8 neighboring cells changes by about 3%. The 
numbers in the plane in the other side are reduced correspondingly. As discussed above, the 
screening charge distribution in the closest plane must be larger and more smoothly varying 
than in this approximation. 

Consider now the effect on the spatial variation of the loop-current order parameter T 
on reducing it to zero by the above effect in the metallic cells close to the muon. The 
amplitude of the order parameter can only recover to its full value over the characteristic 
coherence length £ of the order parameter T(r) beyond the cells where T vanishes. The 
spatial dependence of T is determined by 



where Fq is the free energy of the homogeneous loop-order phase, with the boundary con- 
dition that T is zero in the lattice cells close to the muon. This equation also occurs in the 
context of the variation of the order parameter in the core of a vortex except that here the 
boundary condition T(r=0) = is replaced by T(r) = for r < a , the radius in which the 




(5) 
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local charge density is in the insulating range. A dimensional argument suggests that at low 
temperatures £ ~ a(E e /E c i) 1 ^ 2 , where a is the in-plane lattice constant, E e is a typical elec- 
tronic energy (e.g., the conduction electron bandwidth ~ 1 eV) and E^ is the condensation 
energy of the loop-order phase. Even if we take for the latter a value as high as 0.1 eV, £ is 
about 3a. Near the ordering temperature £(T) diverges. 

Thus the loop-current order parameter is suppressed for a number of lattice constants 
beyond the cells where it is truly zero. Since one or more cells adjacent to the muon have 
zero order parameter, in estimating the magnetic field at the fi + site it is reasonable to 
exclude near-neighbor cells as well as the cell immediately below or above the muon. Then 
the estimated field at the fi + site is consistent with experiment as shown below. 

The calculation is done by setting up an orthorhombic lattice, placing a fi + site and local 
dipole moments in a magnetic structure in the unit cell, and summing over the local moments 
out to a given radius. The muon is placed 1 A below the apical oxygen as suggested by some 
experiments [3], although there seems to be no clear consensus on this location [1]. With no 
change in the loop-current order parameter due to screening of the fi + charge, the magnetic 
field at the /i + site for a magnetic order with a moment of 0.05 \xb per triangular plaquette 
placed at the centroid of the triangles is about 240 G for L^-zSr^CuO^ (a value of 41 G is 
obtained from a more precise calculation of the field due to the loop current pattern [4J.) If 
we now exclude 9 cells in the immediate vicinity of the fi + site in each of the two planes, i.e., 
ao + £ = 2 lattice parameters (7.6 A), the calculated n + magnetic field is about 3.5 G. When 
the next neighboring set of cells are excluded, i.e. ao + £ = 3 lattice parameters (11.4 A) 
the field is about 1 G, while for ao + £ = 4 lattice parameters (15.2 A) it is about 0.45 G. 
Given the uncertainty of the muon site, moving it 15% closer to the apical Oxygen leads to 
a field of 4.35 G, 1.6 G and 0.69 G when screened to two, three and four lattice spacings 
respectively. The above discussion does not take into account the small correlation length 
found for the loop order in La 2 _ x Sr ;c Cu04 [h]], which could further reduce the estimated 
field. 

For the case of YBa2 0^06+5 we find an unscreened /i + field of about 290 G for the 
muon at 1 A below the apical oxygen and the same loop-current magnetic moment as 
above. The loop-current reduction due to screening reduces this to ~3 G when only the 
cells immediately above and below the muon are excluded, and to ~0.5 G when nearest- 
and next-nearest neighbor cells are also excluded. 
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The precise values obtained from these calculations vary with details that are not accu- 
rately known, such as the muon position. The principal point is that the /i + field due to loop 
currents is reduced by more than two orders of magnitude due to the inescapable effects of 
charge screening. 

For metallic doping levels the lattice cells near a fi + site with charge density close to 
that of the insulator will acquire local spin moments at Cu sites. These antiferromag- 
netically-coupled Cu spins will not freeze because of thermal and quantum fluctuations 
in such small regions, but the /i + -induced Cu-spin cluster should act as a paramagnetic 
center and thus might be observable in /iSR experiments. At temperatures well above 
the superconducting transition /^SR measurements give some evidence for weak dynamic 
relaxation in YBa 2 Cu30 6+a: and La 2 _ a; Sr a ,Cu04 2(1 21]; in La 2 _ x Sr x .Cu04 a ~2 G field 
fluctuating at a very slow rate (0.1 MHz) was inferred 21] . Although a Cu-spin cluster might 
be expected to fluctuate slowly, at present its dynamics are not sufficiently well understood 
to allow comparison with experiment. 

Such a cluster should, however, have a paramagnetic response appropriate to its net 
moment, which should lead to a muon Knight shift proportional to 1/T when an external 
magnetic field is applied. This temperature dependence is quite different from that expected 
from bulk susceptibility measurements. Accurate muon Knight shift experiments have not 
yet been carried out [3]. Transverse- field /iSR experiments have, however, detected inhomo- 
geneous static distributions of fi + fields with widths proportional to the applied field 20]. 
E.g., muons see a field-induced distribution of inhomogeneous magnetic fields with a spread 
of about 2 G at a field of 3 T. Suppose there is a 1 [Mb moment produced in the insulating 
cell above and below the muon. With a Curie susceptibility, it produces an extra magnetic 
field at the /z + site of 0(^bH /T)fiB/\r\^ , where |r| is a typical distance to the muon of 
about 5 A and H is the applied field. At a field of 3 T and a temperature 100 K, this 
field is about 100 G. If we allow a distribution of the detailed position of the muons due 
to local strains, a field distribution with a width oc H/T is expected in the normal state. 
Such behavior is consistent with the observations [20]. The observed width of about 2 G at 
H/T = 3 x 10~ 3 T/K is also consistent with our ideas. To test them further, a fi + Knight 
shift oc 1/T should be looked for. 

We are grateful to P. Bourges, A. Kapitulnik, and J. Sonier for useful discussions. This 
work was supported in part by the U.S. NSF, Grant no. 0422674. 
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